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■ Abstract: 

The deconfinement transition in SU (2) gauge theory and the magnetization transition 
Pn| I in the Ising model belong to the same universality class. The critical behaviour of the 
^ I Ising model can be characterized either as spontaneous breaking of the Z2 symmetry of 

^ ' spin states or as percolation of appropriately defined spin clusters. We show that decon- 

O . finement in SU{2) gauge theory can be specified as percolation of Polyakov loop clusters 
^ . with Fortuin-Kasteleyn bond weights, leading to the same (Onsager) critical exponents as 
the conventional order-disorder description based on the Polykov loop expectation value. 
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Colour deconfinement is a well-defined phase transition in finite temperature SU{2) 
gauge theory; the expectation value of the Polyakov loop serves as an order parameter 
Qh! determining the onset of a spontaneous breaking of a global Z2 symmetry |l|, |^. The 

resulting critical behaviour belongs to the universality class of the Ising model, conjec- 
tured on the basis of effective theories ^ and confirmed by lattice studies f^. In full 
QCD, dynamical quarks act as (small) external field; hence the conventional formalism of 
spontaneous symmetry breaking does not define an order parameter, and it is not clear 
5^ I if deconfinement remains a genuine phase transition. It therefore seems helpful to con- 
sider an alternative approach to critical behaviour in the Ising model, which may be more 
readily generalizable to full QCD. 

The magnetization transition in the Ising model (in the absence of an external magnetic 
field) can be described either as the spontaneous breaking of the Z2 symmetry of the 
theory by spin states or as percolation of cluster states appropriately defined in terms 
of the basic spin-spin interaction §]. More specifically, the partition function on a two- 
dimensional lattice of sites is conventionally defined as sum over all possible spin states 
{s} = {s, = ±l Vz = l,...,L2}, 



with the Boltzmann weight 



z{T) = Y.K{s}) (1) 

{s} 

b = exp{{J/T)Y,s,s,}, (2) 
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where T is the temperature and J the spin-spin couphng strength; the sum (i, j) runs 
over next neighbour spin pairs, with i < j. Equivalently, it can be written as sum over 
all 2^ clusters states G, specified in terms of I bonds forming C connected clusters on a 
lattice consisting of E links, 

Z(T) = ^2S'. (3) 

G 

where 

t; = (exp{2(J/T)}-l). (4) 

defines the bond weight. Note that Eqs. (|l|) and (j^) differ by a T-dependent factor related 
to the zero in energy. 

In the spin formulation, the order parameter is given by the spontaneous magnetization 
m(T), defined as the average spin per lattice site. 

Its behaviour near the critical point is governed by the critical exponent f3, with 

m(T) ~ (T, - Tf, T<T,- (6) 
the divergence of the corresponding susceptibility is with 

Xm{T)^{s')-{\s\f^\T-nr. (7) 

determined by the exponent 7. 

In the percolation formulation, the order parameter becomes the percolation strength 
P(T), defined as the probability that a randomly chosen site in the thermodynamic limit 
belongs to an infinite cluster, 

P(T) ~ (T, - T)^ T<T,. (8) 

The corresponding susceptibility is the average cluster size S(T), excluding percolating 
clusters, and its divergence is governed by 

S{T) \T -T,\-^ (9) 

in the vicinity of the critical point. 

Spin and cluster formulations thus provide two equivalent ways to specify the critical 
behaviour of the Ising model. While the spin version is based on the onset of spontaneous 
breaking of the global Z2 symmetry of the Ising Hamiltonian, the cluster version uses the 
onset of percolation of clusters whose Fortuin-Kasteleyn bond weights are also determined 
by Ising dynamics. 

Since for a fixed space dimension d, the critical behaviour of the Ising model and of 
finite temperature SU (2) gauge theory are in the same universality class, one may expect 
that deconfinement can also be formulated as percolation. The aim of this paper is to 
show that for d = 2 and for a specific lattice regularization, this is indeed the case. The 
generalization to ci = 3 appears straight-forward and is in progress. 
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To implement the cluster formulation, we use the droplet approach introduced by 
Coniglio and Klein ||^. For the Ising model, this method generates equilibrium configu- 
rations using the Boltzmann weights and then defines clusters as regions of parallel spins 
connected by bonds, using Fortuin-Kasteleyn bond weights 

p = 1 -exp{-2(J/T)}. (10) 

The use of both Boltzmann and bond weights can be avoided, and there exist implemen- 
tations using the Fortuin-Kasteleyn bond weights only [10| - |jll|]. In our context, the 



Coniglio-Klein implementation seems preferable mainly for computational reasons. 

In the Ising model, clusters are thus defined as regions of parallel spins connected by 
bonds, with the probability for bonding given by Eq. (p^O]). In two-dimensional SU{2) 
lattice gauge theory, the underlying manifold becomes a A^^ x A^^ lattice, with A^^ sites in 
each spatial and Nt sites in the temperature direction. The Ising spins Sj = ±1 at sites 
i = 1, are replaced by Polyakov loops Lj, 

L,~|Trn[/(,^.,.+i)|, (11) 

where U(^i.r,T+i) are SU{2) matrices associated to the link at spatial site i connecting the 
temporal planes r and r + 1 . 

The matrix product in Eq. (|ll|) becomes a loop closed by periodicity in the temperature 
direction. We thus replace the discrete spin values Sj = ±1 by spins of continuous size 
Lj = ±|Lj| at each lattice site i. Clusters are now defined as regions of like-sign Polyakov 
loops (say Lj > 0) connected by bonds distributed according to the bond weight 

Pij = 1 - exp{-2ftLjLj}, (12) 

where k is determined by the dynamics of SU{2) gauge theory; it corresponds to the J/T 
in the Ising weight Eq. (??). 

The identification of the underlying dynamics as SU (2) gauge theory is contained in 
the K in Eq. (0). A general solution for this does not seem to exist so far. In the context 
of universality studies P, ^, the SU{2) action is formally reduced to an effective action 
in terms of the Polyakov loop variables Lj, so it should be possible to identify k in such 
a framework. We shall here make use of strong coupling calculations |T^, in which 



n ^ CN^f (13) 

was found to provide a good approximation (90% accuracy) for Nt = 2; here (3 = {4:/g'^) 
denotes the coupling parameter in the SU{2) lattice action. 

With this, the Polyakov loop percolation problem is fully defined, and we proceed with 
a lattice study of the percolation strength P{l3) and the cluster size S{(3). Our analysis 
is based on four sets of data on A^^ x 2 lattices, with Ns=Q4, 96, 128 and 160. The 
simulations were carried out on workstations for A^^ < 128 and on a Cray T3E (ZAM, 
Jiilich) for A^s=160. The updates consist of one heatbath and two overrelaxation steps. 
For the 64^ x 2 and 96^ x 2 lattices we evaluated configurations every six updates, for 
128^ X 2 and 160^ x 2 every eight updates, measuring in each case -P(/5) and S{/3). The 
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Figure 1: Cluster size as function of (3 near the critical point l3c{L) (dashed line). 
10 



§ 8 

T3 

CD 
CD 

CD 
03 

^ A 
CD 4 

•D 

CD 




P 

S 



3.43 3.435 3.44 3.445 3.45 3.455 3.46 

P 

Figure 2: distributions for scaling fits of -P(/3) and S{j3). 



expectation value of the Polyakov loop gives /3c{L) ~ 3.464 as the critical point on x 2 
lattices in the limit A^s oo fl^, giving us an idea of the range to be studied. 

A first scan for values 3.1 < (3 < 3.5 leads to the behaviour of S{(3) shown in Fig. 1. It 
is seen that S{/3) peaks slightly below Pc{L); with increasing A^^, the peak moves towards 
Pc{L) and the peak height increases. 

In a second step, we carried out high-statistics simulations in a narrower range 3.410 < 
f3 < 3.457 around the transition. In general, we performed between 30000 and 55000 
measurements per /3 value, with the higher number taken in the region of the interval 
closest to the eventual critical point. To obtain the behaviour in the limit Ng ^ oo, we 
applied the density of states method (DSM) [jl4|. This method generates for each lattice 
size by interpolation further values of P{l3) and S{(3). The percolation critical point 
(3c{P) and the critical exponents are then determined through the A^^-dependence of the 
observables. At (3c{P), fits of the form P{(3) ~ N-i^/" or S{(3) ~ N]/" should lead to a 
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minimal |]T5[. Fig. 2 shows the x^/d.f. resulting from fits of logP and logS* vs. log A^^ 
at each value of /5 in the range 3.43 < (3 < 3.46. The two curves show pronounced 
minima with remarkable overlap. 

In Table 1 we show the results for 95% confidence level, comparing the critical Pc{P) 
for Polyakov loop percolation to the f3c{L) from spontaneous symmetry breaking ||13[| , 
and the percolation exponents to the Onsager values for the two-dimensional Ising model 
{(3=1/8, 7=7/4, z/=l). The critical exponents are seen to be in excellent agreement. 
Since z/ is determined by taking derivatives of the interpolated curves of P{P) and S{(3) 
with respect to (3, the precision of its determination is reduced in comparison to that 
oi (3/v and 7/z/. The critical values of j3 are very close, but they do not overlap within 
errors. In view of the approximate solution for used here, small deviations are not 
unexpected. Monte Carlo renormalisation group techniques may well allow a more precise 



determination of the bond weight (p!2D . [|T6 
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Table 1: Critical parameters for SU{2) lattice gauge theory with A^^ = 2. 

We thus conclude that for the specific case considered here, two-dimensional SU{2) 
lattice gauge theory with Nt = 2, deconfinement can indeed be specified through Polyakov 
loop percolation. To make this more general, further work is necessary. The most crucial 
open problem is certainly the general form of k in Eq. {\L2^ . In particular, it is not 
clear why the general SU{2) action should lead to a nearest-neighbour form, and how 
the 'temperature' k can be defined more generally, away from the strong coupling limit 
Nt = 2. For the latter problem, the use of Monte Carlo renormalisation group techniques 
p!7| may be of help. If and when these problems are solved, the corresponding percolation 



study has to be carried out for three space dimensions and for SU{^) gauge theory as 
well, where the transition becomes first order. 

Finally we return to the motivation mentioned at the beginning. For spin systems, 
percolation in coordinate space remains as critical phenomenon even in the presence of 
an external field although this criticality is now not connected any more to singu- 
lar behaviour of the partition function. In principle, this could provide the basis for a 
deconfinement order parameter in full QCD with dynamical quarks [|19] . 
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